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Abstract 

This paper examines the Hausdorff dimension of the level sets f^ 1 (y) of 
continuous functions of the form 

oo 

f(x) = 2- n u n {x)^{2 n x), < x < 1, 

n=0 

where 4>(x) is the distance from x to the nearest integer, and for each n, uj n is 
a { — 1, l}-valued function which is constant on each interval [j/2 n , (j + l)/2 n ), 
j = 0, 1, . . . , 2 n — 1. This class of functions includes Takagi's continuous but 
nowhere differentiable function. It is shown that the largest possible Hausdorff 
dimension of f (y) is log ((9 + \/105)/2) / log 16 « .8166, but in case each ui n 
is constant, the largest possible dimension is 1/2. These results are extended 
to the intersection of the graph of / with lines of arbitrary integer slope. 
Furthermore, two natural models of choosing the signs <jj n [x) at random are 
considered, and almost-sure results are obtained for the Hausdorff dimension 
of the zero set and the set of maximum points of /. The paper ends with a 
list of open problems. 
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1 Introduction 



Takagi's continuous nowhere different iable function, shown in Figure Q^a) 
by 



is defined 



where <fi{x) = dist(x, Z), the distance from x to the nearest integer. A great deal has 
been written about this function since its initial discovery in [23]; see recent surveys 
by Allaart and Kawamura [6] and Lagarias [15] for an overview of the literature. In 
the past few years, interest has focused mainly on the level sets L(y) := {x G [0, 1] : 
T(x) = y}, which have been shown to possess a rich structure. For instance, L(y) is 
finite for Lebesgue-almost every y [11] , and can have any even positive integer as its 
cardinality (3] . However, the "typical" level set of T (in the sense of Baire category) 
is uncountably large [I]. These uncountable level sets can be further differentiated 
according to their Hausdorff dimension. Kahane [32] showed that max x T(x) = 
2/3, and L(2/3) is a Cantor set of Hausdorff dimension 1/2. De Amo et al. [7] 
recently proved that 1/2 is the maximal Hausdorff dimension of any level set, settling 
a conjecture of Maddock [18J, who had earlier obtained a bound of 0.668. Two 
interesting papers by Lagarias and Maddock [161 ELY] use novel notions of 'local level 
sets' and a 'Takagi singular function' to establish several properties of the level sets 
of T. For instance, it is shown in [T7] that the set of ^/-values for which L(y) has 
strictly positive Hausdorff dimension is a set of full Hausdorff dimension 1. 

In this article we examine the level sets of a class of generalized Takagi functions 
of the form 

OO , s 

/(*)=£ ^0(2"*). 0<x<l, (1) 

71=0 

where 



-M}, 



constant on each interval 



J_ 

2 n> 
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0,1, 



so that the terms of the series are continuous and hence, / is a continuous function. 
We denote by % the class of all functions / of the form ([I]), and by T c the subclass 
of those / in % for which u n is constant for each n. The class T v was investigated 
in detail by Abbott, Anderson and Pitt [T], who denoted it by K* dl because of its 
relationship with Zygmund's class A* of quasi-smooth functions. But whereas [T] 
studies the class % from the perspective of abscissa or x-values, our focus here is on 
ordinate or y- values. 

Several members of T v have featured in the literature. These include the alter- 
nating Takagi function (e.g. [TJ [13]), for which u n = (— l) n and which hence lies in 
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Figure 1: The Takagi function (left) and the Gray Takagi function 



%; and the Gray Takagi function [H], for which u n (x) = r n (x), where r n denotes 
the nth Rademacher function defined by 

r n (*):=(-l) L2 "* J , n = 0,l,2,... (2) 

The Gray Takagi function is shown in Figure []Jb). Another example is the function 
T 3 of Kawamura [TO], which has uj n (x) = r\(x) ■ ■ ■ r n (x). All members of T v are 
nowhere differentiable; Billingsley's proof [TO] for the Takagi function extends easily. 
Furthermore, their graphs have Hausdorff dimension one [8J. All functions in the 
subclass T c are symmetric with respect to x = 1/2, and their level set structure is 
similar to that of T: Lebesgue-almost all level sets are finite, but the "typical" level 
set is uncountably infinite; see [5]. Whether these properties hold for the wider class 
T v remains unsolved. 

This paper concerns the Hausdorff dimension of the level sets 

L f (y) := {xe [0,1] :f(x)=y}, fe%, y G R. 

For a Borel set E C R, we denote the Hausdorff dimension of E by dimn E. 

1.1 Sharp upper bounds 

The first half of the paper gives sharp upper bounds for diniff Lf(y), both for / e % 
and for / G %. 

Theorem 1. For any f G T c , we have 

maxdiniff Lf(y) = 1/2. (3) 
y 

The Gray Takagi function satisfies even though it does not lie in T c . 
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Theorem 2. Let f be the Gray Takagi function, 



/(*) = £ ^#2"*). W 

n=0 

where r n is the nth Rademacher function. Then max y dim// Lf(y) = 1/2. 

For general f £ 71, however, the Hausdorff dimension of the level sets of / can 
be much greater: 

Theorem 3. Let a := (9 + vT05)/2, and put d* v := log a/ log 16. We have 

max max dim# = d* « .8166. (5) 

By embedding the graph of / affinely into the graph of a suitable function g G 7~ v , 
we can extend the above results to the intersection of the graph of / with arbitrary 
lines of integer slope. 

Corollary 4. For m,b £l, let £ m ^ denote the line with equation y = mx + b. 
(i) For each f G% and each m G Z, 

maxdim//(Graph(/) f\C m ,b) = 1/2. 

(m) Let d* be as in Theorem^ For each m G Z, 

maxmaxdiHifl-fGraphff) n£ m b) = rf,*,- 

We prove Theorems [U [2] and [3] and Corollary |4] in Section [3], by modifying and 
extending the method of De Amo et al. [7] . The idea is to consider the intersection 
of the graph of the partial sum f 2n of / (defined in ([H]) below) with a suitably chosen 
horizontal strip [0, 1] x J n , where { J n } is a nested sequence of intervals shrinking to 
{y}. We then derive a system of linear recursions for the number of line segments 
in these strips (differentiated according to their slopes). An added complication is, 
that the coefficients in these recursions are dependent on n. Thus, for instance, the 
greater part of the proof of (jSJ) consists in determining the joint spectral radius of a 
certain pair of 3 x 3 nonnegative matrices. It appears to be a lucky coincidence that 
this can be done exactly. 
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1.2 The random case: zero sets and maximal sets 

Perhaps most interesting is the case when the signs u n are chosen at random. We 
consider the two most natural schemes here. Let Z + denote the set of nonnegative 
integers. 

Model 1. (Random choice from T c ) Assume l)q,l)x,... are independent and 
identically distributed (i.i.d.) with P(u n = 1) = p = 1 — P(u n = —1), where 
< p < 1, and set co n (x) := u n . 

Model 2. (Random choice from T v ) Assume {oo n ,j '■ n G = 0, 1, . . . , 2 n — 1} 
are i.i.d. with P(co n ,j — 1) — p — 1 — P(uj n j = — 1), and set 0J n (x) := u n j if 

x e b'/2 n , (j + l)/2 n ). 

In either model, set q := 1 — p. Since it seems difficult to treat the level sets in 
full generality, we focus on two special cases: the zero set and the set of maximum 
points. Let 

do ;= M(l + V5)/2) m 
log 4 

For the zero set of /, we have the following results. 
Theorem 5. Assume Model 1. 

(i) The zero set Lf(0) is finite with probability 

P(#L/(0) < oo) = 1 - min{p/q, q/p}, 
and given that Lf(0) is infinite, dim// Lf(0) > a.s., for each p G (0, 1). 

(ii) If p = 1/2, then dim H Lf(0) < d a.s. 

Theorem 6. Assume Model 2. If p = 1/2, then dim#L/(0) = do a.s. 

Which functions / G %, specifically, satisfy dim H L/(0) = d 7 Again the Gray 
Takagi function provides an example. 

Proposition 7. Let f be the Gray Takagi function. Then dim# Lf(0) = do. 
Next, for / G %, define 

Mt := max f(x), 

1 x G [0,l] 

and 

M f := {x G [0, 1] : f(x) = M f } = L f (M f ). 

Note that / <T for all / G %, so Mf < 2/3. The following theorem was proved in 
[2] and is included here for comparison with Model 2. 
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Theorem 8. Assume Model 1. 



(i) If p > 1/2, then the distribution of Mf is singular continuous and dmif/.M/ = 
1 - (2P)- 1 a.s. 

(ii) If p < 1/2, then the distribution of Mf is discrete and M.f is finite a.s. 

In fact, the paper [2] specifies the distributions of Mf and the cardinality of M.f 
in considerable detail under the assumption of Model 1. The analysis appears to be 
much harder for Model 2, and we describe here only what happens when p > l/y/2. 
Note the contrast with the previous theorem. 

Theorem 9. Assume Model 2. 

(i) The probability that Mf attains the maximum possible value of 2/3 is 

P(M f = 2/3) = max( 2p ~ 1 ,0 

I V A 

(ii) Ifp > l/y/2, then 

log(2p 2 ) 

duRaMf = ^T a - s - } 

and the distribution of Mf is discrete, supported on the set {^YljeA^ • ^ ^ 
Z+, #(Z+\A) <oo}. 

It seems plausible (by monotonicity considerations) that dim// Aif = a.s. for 
p < 1/ \/2, but the author has not been able to prove this. 

The results for the random case are proved in Section HI Theorem O is proved 
by casting the zero set as the attractor of a Mauldin- Williams random recursive 
construction and by using properties of hitting times in a symmetric simple ran- 
dom walk. The proof of Theorem [5] requires a different approach. For the upper 
bound, we use the Perron- Frobenius theorem to construct a sequence of positive su- 
permartingales which lead, via the Martingale Convergence Theorem, to successive 
upper bounds for dim#L/(0); we then prove that these upper bounds converge to 
do. For the lower bound we identify a particular combination of line segments, called 
a "Z-shape" , which will appear around the x-axis in the step-by-step construction of 
the graph of / with probability one given that Lf(0) is nonfinite. We show, via the 
law of large numbers, that the number of Z-shapes grows exponentially fast almost 
surely once a Z-shape appears. This, along with some additional observations, gives 
the lower bound. Finally, the basic theory of branching processes is used to prove 
Theorem |9j 

There are many natural questions still unanswered; Section |5] lists some of them. 
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2 Preliminaries 

The following notation will be used throughout. For an interval J C [0, 1], \ J\ denotes 
the diameter of J, whereas the cardinality of a discrete set A is denoted by #A. For 
n G Z + and / defined by (OQ), put 

71-1 



fc=0 



Then /„ is piecewise linear, and the right-hand derivative of f n at any point x G [0, 1) 
is 

71-1 

fn( x ) = ^2uJ k (x)r k+1 (x), 

k=0 

where r n is defined as in It is clear that — fn( x )\ — 1 f° r an x an d all 

n, and so fn +2 ( x ) ~ fn( x ) e { — 2, 0, 2}. In particular, f^x) is always even. 

Another important observation is that f(k/2 n ) = f n (k/2 n ) for integer k and 
n G Z + . We need two more elementary facts about the functions f n : 

Lemma 10. For each n G Z +; \ f — f n \ < 2~ n . 

Proof. By the periodicity of 0, 



fc=n 



2 _ n ^ 2 - m 0(2 m 2 n x) = 2~ n T(2 n x mod l) . 



m=0 

Thus, the lemma follows from the bound T(x) < 2/3. □ 
Lemma 11. For n G Z + and integer j , 4 n /2n(j/2 2n ) is an even integer. 
Proof. We can write 

4n /- (Jz) = E ^ ( 2^) 2 "<^ (2^) • ( ? ) 

Now for each integer m > 1, 2 m 0(j/2 m ) is an integer with the same parity as j. To 
see this, note that <j)(j/2 m ) is either j' /2 m or (2 m — j')/2 m , where j' = j mod2 m , and 
j, j' and 2 m — j' all have the same parity. Thus, the terms in the right hand side of 
(J7J) are either all even or all odd. Since there is an even number of them, their sum 
is even. □ 
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We now introduce the closed intervals 
j j + 1 



In ' On 



nGZ + , j = 0,l,...,2"-l, 



J, 



n,k 



2k 2k + 2 



n G Z+, fceZ. 



4™ 7 4« 

Let J° ■ denote the interior of I n j. 

3 Universal upper bounds 

Following De Amo et al. [7], we divide the strip [0, 1] x R into closed rectangles 
R n ,j,k ■= hnj x J„, fc , n G Z+, < j < 2 2n , k G Z. 



Let 



Put 
and 



T'n.fe := {j : -Rn,j,fc intersects the graph of /}, 
Kk ■= {j:f2n = 2k/A n on I 2nd }, 
E\ k := {j : the interior of -R n j,fc intersects the graph of /W}- 

:= #-E n , fe , iV° fc := #£° fc , iV^. := #£ n , fc , 



M n := maxiVnjfc, 

k&L 



M° := maxiV; 



n,ki 



Ml := maxN^ k . 



It is important to note that if j G £?„ fc , then R n ,j,k intersects the graph of f2 n in a 
nonhorizontal line segment, in view of Lemma ITT1 Another key observation is that 
if R n ,j,k intersects the graph of /, then the graph of jm must intersect R n ,j,k-i U 
Rn,j,h U Rn,j,k+i- This follows easily from Lemma fTUl Hence, for each k, 



N n ,k < N^ k + A° >fe+1 + iV^ jfc _ 1 + N n>k + N n>k+1 , 



so that 



M n <2M° + 3M*. 
Proof of TheoremUl Let / G Tc- We claim that 



and 



M n ° +1 <max{2M n °,M^}, 



Mi +1 <2M n ° + M^. 



(8) 

(9) 
(10) 
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To verify (jUJ) and (ITUj) . we consider four cases regarding co> 2n and ci^n+i- Let s n j be 
the slope of /2 n on l2njj and put s n+1 j := (s n+ x,4j; s n+i,4j+i) ^+1,4,7+2, ^+1,47+3)- 

Case 1: uj^n = <^2n+i = 1- If j £ then s nj - = and hence s„ +lj - = 

(2,0,0, -2). Thus, 4j + 1 and 4j + 2 lie in £° +14fc+1 , while 4j and 4j + 3 lie in 
+i,4fe- ^ .? ^n,fe' then /2n+2 is monotone on /2n,j- 111 particular, if s n j — 2, 
then s n+ ij = (4, 2, 2, 0) and hence, 4j + 3 G E® +1 4k+4 . Similarly, if s n j = —2, then 
4j G E® +1 Ak+4 . If \s n j\ > 4, then no component of s n+1 j is equal to 0. We see that: 

• For k' — 4k + 1, each j in fc contributes two members to E® +1 k , while each 
j in k contributes none. Hence 

K + iAk + i = 2< fc . (11) 



For fc' = 4fc, each j in fe or j in E\ k contributes no members to E° +1 k , while 
each j in E\ k _ x contributes at most one. Hence, 

K +1Ak < N^ k _ v (12) 



For k' = 4k + 2 or k' = 4k + 3, neither j in fc nor j in fe contributes to 



E n+ i t y, so 



-^n+l,4fc+2 ~~ ^n+lAk+3 ~~ ^. (13) 



From dnj, flT2) and jSJ), © follows in Case 1. 

The verfication of fflUj) is simpler: It is clear that each j £ E® k contributes at 
most two members to any set E* +1 k , and each j G E\ k contributes at most one by 
the monotonicity of /2n+2 on hnj- From this, (jTUJ) follows. 

Case 2: Co> 2 n = 1, ^2«+i = — 1- Now if j E E® k we have s n+1 j = (0, 2, —2, 0) and 
so 4j and 4j + 3 lie in £° +Mfc , while 4j + 1 and 4j + 2 belong to E l n+lAk . If j G E l n k 
then /2n+2 is monotone on hnj, so j contributes at most one member to any set 
En + i t k>- ^ iri particular s n ,j = 2, then s n+1 j = (2,4,0,2) and so 4j + 2 G -E° + i i4fc+3 - 
Similarly, if s n j = —2, then 4 j + 1 G E^ +lik+z . If |s nj | > 4, then s„ +lj - does not 
have any zero components. We conclude that: 

• For k! = 4k, each j in E® k contributes two members to E° +i k , while j in E\ k 
or j in E\ k _ x contributes none. Hence 

N° = 2iV n %. (14) 
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For k! = 4k + 3, each j in E® k contributes no members to E® +1 k , and each j 



in E\ k contributes at most one. Hence 



iV n ° +li4fc+3 < < fc - (15) 

Finally, for k! = 4k + 1 or k! = 4k + 2, neither j in E^ k nor j in E\ k contributes 
to E° n+l k ,, so 

^n+lAk+l = ^n+l,4fc+2 = 0- (16) 



From ( fl4l) . ( jl5l) and ( 1161) . (Ej) follows in Case 2. And ( TTOT) follows in the same way as 
in Case 1. 

Case 3: uj2n = — l)^2n+i = 1- This case is similar to Case 2, by symmetry 

Case 4 : ^2n — ^2n+i — — 1- This case is similar to Case 1, by symmetry. 

Using (|9]), ( TTTJT) and the initial conditions Mq = 1 and Mq = 0, it is straightfor- 
ward to verify inductively that 

M° < 2 n , Ml < 2(2 n - 1). 

Hence, by ©, M n < 2 n+1 + 6(2 n - 1) < 2 n+3 . It follows that for each y, L f (y) is 
covered by at most 2 n+3 intervals of length 4~ n . Thus, 

(n + 3)log2 1 



dim// L f (y) < dim B L f (y) < lim 



n log 4 2 



It remains to show that for each f E T c there is an ordinate y with dim# Lf(y) > 
1/2. The equalities in (TTTj) and (fT4l) show (by induction) that there is a sequence of 
ordinates y n := k n /4 n , n G Z+, such that AT° fc = 2 n for each n. Specifically, k = 0, 
and 

{4fc n + 1, in Case 1, 
4k n , in Case 2 or 3, 
4k n — 1, in Case 4. 

Let := U{-^2n,j : j G -^nfc }■ The sets {K n } are compact and nested, so the 
intersection K := f]^ =Q K n * s a nonempty compact set. The intersection of K n+ \ 
with any interval I := l2 n j in consists of two (non-overlapping) sub intervals of 
one-fourth the length of /. Hence, by the basic theory of Moran fractals, dim# K = 
log2/log4 = 1/2. Now let y := lim^oo y n , which exists since \y n +i — y n \ < 4~( n+1 \ 
If x G K, then f2n{%) = Vn f° r eac h n, and taking limits gives f(x) = y. Hence, 
Lf(y) D K, and so dim#L/(?/) > 1/2. □ 
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Proof of Theorem [H Let / be the Gray Takagi function, and let s n ^k be the slope 
of /271 on I 2n ,k- The Gray Takagi function has the special property that for each n, 
s n fi = In and |s n ,fc+i — Sn,fc| = 2 for all k, as is easily checked by induction. This 
implies that s n j = 2n + 2j (mod 4) for all n and all j. Thus, we have the following 
substitutions 

n even: ^(2,0,-2,0) n odd: ^(0,-2,0,2) 

2 ^ (2,0,2,4) 2 ^(4,2,0,2) 

-2 ^(-2,-4,-2,0) -2 ^(0,-2,-4,-2). 

Now we can infer that each horizontal line segment in the graph of f2 n spawns two 
horizontal line segments at different levels in the graph of f 2n +2, and each line segment 
of slope 2 in the strip [0, 1] x J n ^ spawns one horizontal line segment at a different 
level than that spawned by a line segment of slope —2 in the same strip. Of course, 
if \ s n,j\ > 4, then f 2n+2 is strictly monotone on I 2n j- Since f 2n (fy = /2n(l) = 0, 
exactly half the k line segments intersecting the interior of the strip [0, 1] x J n k 
have positive slope. From these observations, we conclude that 

M° +1 < Ml + l -Ml M x n+1 < 2M n ° + Ml 

With the initial values M ° = 1 and M X = 0, it follows by induction that M° < 2 n_1 
and M^ < 2 n for all n > 1. As in the proof of Theorem [TJ this implies that 
dim h Lf(y) < 1/2 for all y. 

To complete the proof, we will show that dim^ Lf (2/5) = 1/2. Set y\ := 1/2 and 
K\ := [0, 1/2]. The left half of the graph of f'2 consists of two line segments with 
slopes (2, 0). Denote the union of these two line segments by X. In the graph of / 4 , 
this figure is replaced by 8 line segments with slopes (4, 2, 0, 2, 0, —2, 0, 2). See Figure 
[2j Observe that the third and fourth of these 8 line segments together form a similar 
copy of X, reduced by 1/4 and rotated by 180°; and similarly for the seventh and 
eighth. Moreover, these two copies of X sit at the same height in the graph of f\, with 
their horizontal parts at height y 2 := 3/8. We put K 2 := [1/8, 1/4] U [3/8, 1/2]. Next, 
in the transition from f± to fa, each of these two rotated copies of X, now having 
slopes (0, 2), is replaced by a piecewise linear figure with slopes (2, 0, —2, 0, 2, 0, 2, 4). 
This new figure contains two smaller copies of the original figure X, rotated back 
and sitting at equal height in the graph of fa; see the second part of Figure EJ Thus, 
overall, the graph of /g contains (at least) four similar copies of X with the same 
orientation, and with their horizontal part at height 2/3 := 3/8 + 1/32. Let K 3 be the 
(closure of) the union of the projections of these four copies of X onto the x-axis. 

It is not hard to see that this pattern continues: Each copy of X in the graph 
of f 2n induces two similar copies of X in the graph of f 2 n+2, reduced by 1/4 and 
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— > zoom m 



\ 



/4n-2 



4n 



Figure 2: The substitution (2, 0) h-> (4, 2, 0, 2, 0, -2, 0, 2) when going from / 4n _ 2 to 
/ 4n , and the substitution (0, 2) h-> (2, 0, —2, 0, 2, 0, 2, 4) when going from / 4n to /4 n +2, 
in the construction of the Gray Takagi function. 

rotated 180°. Both copies sit at the same height in the graph of /271+2, with their 
horizontal part at level y n+ i = (1/2) XT=o( — V^) 1 - Thus, by induction, the graph of 
/2n+2 contains 2 n copies of X, equally oriented with their horizontal parts at level 
Un+i- Let K n+ i be the closure of the union of their projections onto the x-axis. 

The sets {K n } are nested; put K := {~)™ =1 K n . Note that y n — > 2/5. Now x G K 
implies that for each n, \f2n(x) — y n \ < (1/2)(1/4)™ _1 , and letting n — > oo we get 
f(x) = 2/5. Thus, K C Lf(2/5). Finally, K is a two-part generalized Cantor set 
with constant reduction ratio 1/4, and hence, dim. H K = 1/2. □ 

3.1 Flexible case: construction of an extremal function 

We construct a function / e % and an ordinate y such that dim H Lf(y) > d* v . We 
start with f = 0, and describe inductively how to construct /2n+2 from f 2n for each 
n. At the same time, we construct a sequence of ordinates {y n }, called baselines, 
which converge to the desired ordinate y. We also construct a nested sequence {K n } 
of compact subsets of [0, 1] such that everywhere on K n , the slope of ji n takes values 
in {—2,0,2} only. This sequence {K n } will converge to a compact limit set K, 
which will be a subset of Lf(y) and will have Hausdorff dimension d*. The idea 
of the construction is to build the successive approximants jin i n such a way as to 
maximize the asymptotic growth rate of the number of intervals l2 n ,j contained in 
the level set L hn (y n ). 

Let u n ,k be as in Model 2 in the random case; that is, u n (x) = u nt k if x G I° k . 
Define the functions 



4>n,k( X ) 




if x G I ntk 
otherwise. 
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Set /o = 0, yo — 0, and K = [0, 1]. For n G Z+, assume that /Wi, 2/n and K n have 
already been constructed in such a way that the slope of fin on K n takes values in 
{—2, 0, 2} only. Fix j e {0,1,..., 2 2n — 1}, and note that on hnj, fin+i is obtained 
from f 2n by adding 

^2n,j4>2n,j + k'2n+l,2.j'02n+l,2j + ^2n+l,2j'+102n+l,2j+l- 

Thus, to describe the inductive construction of /2n+2 from f 2n on J 2n) j it suffices to 
describe how we choose the vector of signs w := (uj2n,j,^2n+i,2j,^2n+i,2j+i)- If hn,j 
does not lie in K n , it does not matter how we choose w, but for definiteness, we put 
w = (1, 1, 1). Assume now that l 2n ,j C K n , and consider four cases. Let s denote 
the slope of f 2n on I 2n ,j, so that s G {—2, 0, 2}. 

Case 1: n = (mod 4). In this case we move the baseline up by setting y n +i = 
Vn + ■ If s = 0, we take w = (1, 1, 1); if s — 2, we take w = (—1, 1, —1); and if 
s = —2, we take w = (—1, —1, 1). 

Case 2: n = 1 (mod 4). In this case we do not move the baseline: y n +i = y n . If 
s = 0, put w = (-1, 1, 1). If s = ±2 and f 2n < y n on I 2nJ , we put w = (1, 1, 1); 
otherwise, we put w = (—1, —1, —1). 

Case 3: n = 2 (mod 4). In this case we move the baseline down by setting y n +i = 
y n — |(|) n - If s = 0, we take w = (—1, —1, —1); if s = 2, we take w = (1, —1, 1); 
and if s = —2, we take w = (1, 1, —1). 

Case 4-' n = 3 (mod 4). In this case we do not move the baseline: y n+ \ = y n . If 
s = 0, put w = (1, —1, —1). If s = ±2, we choose w the same way as in Case 2. 

The transition from f 2n to f 2n+2 is illustrated in Figure [3] for Cases 1 and 2. 
Substitutions not shown in the figure are obtained by symmetry. Cases 3 and 4 are 
mirror images of Cases 1 and 2, respectively. 

Now that /2n+2 and y n +\ have been constructed, we define K n+ \ as follows. For 
j = 0, 1, . . . , 2 2n+2 — 1, put the interval hn+2,j in K n+ \ if and only if I 2n+2 j C K n 
and /2n+2(a^) = y n +i for some x G I 2n +2,j- Then on each interval hn+2,j in K n+ i, the 
slope of f 2n+2 lies in {-2, 0, 2}. 

This completes the construction of {/ 2 „}, {y n } and {K n }. Now let y := lim^oo y n , 
K := n^=o-^ n an d / := nm »woo fn- Then / G %, and it is not difficult to calculate 
that y = 8/17. We claim that K C Lf(y). To see this, let x G K and n G N. 
The slope of / 2 n at x (if defined) lies in {—2,0,2}, and there is a point x n with 

-x\< 2~ 2n and f 2n (xn) = Vn- Hence, |/ 2 n(^) - 3/n| < 2 l ~ 2n . Letting n-^oowe 
obtain f(x) = y. 



13 




Figure 3: The substitutions of Case 1 (top) and Case 2 (bottom). The dotted line 
indicates the baseline. Numbers are absolute values of slopes. 

Proposition 12. We have dim H K = d* v , and consequently, dim H Lf(y) > d* v . 

Proof. We shall represent K as a multi-type Moran set. Recall that K n consists of 
those intervals I 2n j on which the graph of f 2n either coincides with the nth baseline or 
else is a line segment with slope ±2 which has an endpoint on the baseline. Let £ be a 
line segment constituting the graph of f 2n above an interval l2 n ,j in K n , and let s be 
its slope. We classify £ to be of type 1 if s = (so that £ lies on the baseline); of type 2 
if s — ± 2 and £ points from the baseline in the direction of the next baseline change; 
and of type 3 if s = ± 2 and £ points from the baseline in the direction opposite the 
next baseline change. More precisely, £ is of type 2 if sgn(/ 2ri — y n ) = sga(y n+2 — y n ) 
on I 2n ,j, and of type 3 if sgn(/ 2 n - y n ) = -sgn(y n+2 - y n ) on I 2n ,j- Thus, the line 
segments of slope ± 2 above the baseline are of type 2 in Cases 1 and 4 above, and 
of type 3 in Cases 2 and 3. For line segments below the baseline, the classification 
is reversed. Say an interval I 2 nj in K n is of type i if the line segment forming the 
graph of f 2n over I 2n j is of type i, for i — 1, 2, 3. 

The transition from K n to K n+ i is summarized by a 3 x 3 matrix = [ruff], 
where denotes the number of type i intervals in K n+ \ contained in a type j 
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A. .A. ,A, A ., A, .A, A ..A..A, 



Figure 4: Graph of a function / G T v attaining the maximal Hausdorff dimension 
d*. The horizontal line represents the level y = 8/17 for which dim# Lf(y) = d* v . 



interval in K n . Letting A and B be the matrices 



A :-- 



2 10 
2 10 
1 



B :-- 



2 1 1 
2 10 
1 



we see from the four cases above that M^ n ' = A if n is even, and M^ n > = B if n is 
odd. This suggests combining the stages two at a time and looking at the transition 
from Km to K2 n +2, which is governed for each n by the product matrix 



M := BA 



6 4 
6 3 
1 



At this point we can eliminate type 3, because intervals of type 3 in Km do not 
intersect Km+2 and hence do not intersect K. Thus, we may delete the last row and 
column of M, and conclude that K is a two-type Moran set with construction matrix 



M 



6 4 
6 3 



The spectral radius of M is a = (9 + V105) /2. Since each basic interval in K 2n +2 is 
1/16 the length of a basic interval in K 2n , the common contraction ratio is 1/16. It 
follows from Marion [19] that dim# K = log aj log 16 = d* v . □ 



The function / constructed above is shown in Figure HI along with the horizontal 
line y = 8/17 which intersects the graph of / in a set of dimension ci*. 
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3.2 Flexible case: proof of the upper bound 

To prove that d* is a universal upper bound for dini^ Lf(y), we again consider the 
numbers N Hjk , iV° k and k defined at the beginning of this section. In this setting, 
however, it is no longer adequate to examine the intersection of the graph of fin with 
a single strip [0, 1] x J n k . Instead, we consider simultaneously two adjacent intervals 
J n ,k-i and J n k , as outlined below. 

Let y G K. be fixed. For each n, y belongs to some J n j, and we may assume that y 
is not the right endpoint of J n j. If 2j/4 n <y< (2j + l)/4 n , put k n = j. Otherwise, 
if (2j + l)/4 ra < y < (2j + 2)/4 n , put k n = j + 1. Now define the numbers 

c« := < fc „, In ■= < fcn _i, u n := Nl kn (n G Z+). 

The dynamics of the triple (c n ,l n ,u n ) depend in a complicated way on y and /. To 
simplify the analysis, we introduce the new variables 

o~ n '■— l n + u n , Tn n '.= max{/ n , "U n }, 

and consider the dynamics of the vector x„ := (c n , <r n , m n )*. For matrices A = [a^-] 
and B = [bij] of equal size, write A < B if a^- < 6^ for all (i, j). Let and F be the 
matrices 





"2 





1 




"2 


1 


0" 


E := 


2 





2 


F:= 


2 


1 







2 





1 




2 





1 



Lemma 13. For each n, either x„ +1 < _Ex n or x n+1 < Fx„. 

Proof. The choice of /c n implies that k n+ i G {4/c„ — 2, . . . ,4k n + 2}. By symmetry, 
the argument for the case k n+ i = 4k n — i is the same as that for k n+ i = 4k n + i 
(i = 1,2), so we have essentially three different cases to consider. Put y n := 2k n /4 n , 
so y n is the center of the combined interval J n ,k„-i U J n ,k„- 

Case 1: k n+ i = 4k n . In this case, y n +i = y n . Each horizontal line segment in 
the graph of f2 n at height y n contributes at most two horizontal line segments at 
height y n +i in the graph of f2n+2, and each nonhorizontal line segment (whether 
above or below height y n ) contributes at most one such horizontal line segment. 
Thus, c n+ \ < 2c n + cr n . Further, each horizontal line segment in the graph of / 2 „ at 
height y n spawns two nonhorizontal line segments (of slopes 2 and —2) either above 
or below height y n +i, while each nonhorizontal line segment generates exactly one 
such segment. Thus, cr n+ i < 2c n + a n . Finally, u n+ i < 2c n + u n and l n+ i < 2c n + l n , 
so m n+ i < 2c n + m n . It follows that x n+i < Fx n . 

Case 2: k n+ \ = Ak n + 1. In this case, y n+ \ = y n + \{\) n ■ Each horizontal 
line segment in the graph of jin at height y n contributes at most two horizontal line 
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segments to the graph of /2n+2 at height y n +i, while each non-horizontal line segment 
contributes at most one. Nonhorizontal line segments which lie entirely below y n 
contribute none. Thus, c n+ i < 2c n + u n < 2c n + m n . On the other hand, a horizontal 
line segment at height y n in the graph of f2 n contributes at most two nonhorizontal 
line segments (of slope 2 and —2, respectively) to the graph of /2n+2 between heights 
y n and y n +i, and none above height y n +i- Finally, each nonhorizontal line segment 
in the graph of j^n that lies at least partially above y n contributes at most one 
nonhorizontal line segment to the graph of /2n+2 in the strip [0, 1] x J n+ i i / Cn+1 , and at 
most one in the strip [0, 1] x J n +i,k n+1 -i- Hence, l n+ \ < 2c n + u n and u n+ \ < u n . It 
follows that a Tl+ \ < 2c n + 2u n < 2c n + 2m n and m n+ i < 2c n + m n . Thus, x n+ i < Ex n . 

Case 3: k n+ i = Ak n + 2. This is the simplest case, because here c n+ \ = 0, 
ln+i < u n and u n+ i < u n . Hence, a n+ i < 2u n < 2m n and m n+1 < u n < m n . It 
follows that x n+1 < Ex n . □ 

As a side note, we can observe from Case 3 above that if y is such that k n+ \ = 
4k n ± 2 for all but finitely many n, then Lf(y) is finite. This will be the case if 
y/2 has quarternary expansion ending in all 2's, i.e. if y = ±(k + l/3)4 _m for some 
k G Z + and m G Z + . In particular, if y = ±1/3, then L/(y) is either empty or 
consists of exactly two points, for each / G %. 

It transpires from Lemma fT3l that we must compute the joint spectral radius of 
E and F. For a set £ of d x d matrices, the joint spectral radius p(S) is defined by 

p(E) := lim sup sup{ || T{T 2 ■ ■ ■ T n \\ 1/n : T, L G E for all i}. 

The value of p(E) is independent of the choice of matrix norm ||.||; see Rota and 
Strang [22]. In general the joint spectral radius is difficult to compute exactly, even 
for sets of just two matrices. While much work has been done for 2x2 matrices (see 
Mossner [21] and the references therein), there are few known examples for larger 
matrices. For our set E = {E, F}, however, it is possible to show through a sequence 
of steps that p(E) = ^fa, and this will give the upper bound in Theorem [3j 

Proposition 14. The joint spectral radius of E and F is given by 

p({E,F}) = ^(FE) = ^, 
where for any square matrix M, p(M) denotes the spectral radius of M. 

Proof. We use the matrix norm ||T|| = l%l> where T = [i^-]. If T has nonnega- 
tive entries, we have the representation ||T|| = 1T1*, where 1 = (1, 1, 1). Denote by 
Vk the set of all products Tx ■ • • T\. where % G {E, F} for i = 1, . . . , k, let Vq := {/}, 
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and put V := Ufclo^- ^ useful first observation is that for M G T 7 , we have the 
implications 



1M 
Ml 1 



{Wi,W 2 ,W 3 ) 

as is easily checked by induction. Since 



E 



6 3 
8 4 
6 3 



«i > max{f 2 ,v 3 }, 
min{wi,W2} > W3, 



6 2 1 
/:;F= 8 2 2 
6 2 1 



(17) 
(18) 



we see by flE} that for each Me?, £ 2 M1* < EFMl*. Thus, if £ 2 occurs anywhere 
in the product T1T2 ■ • -T n , we can replace it with EF without decreasing the norm 
of the product. Hence it suffices to consider products of the form 



F h EF i2 E . . . F jk E ^ k e Nj 3i > 1? 
Lemma 15. Lei S be a matrix of the form 



(19) 



5 



a 
a 
0-6 



a > 6 > 0. 



Then for allM 1 ,M 2 E V, 



lM l SM 2 l t > 0. 



Proof. Let v = (t>i, t>2, ^3) = lMi, and w = (wi,W2,wsY = M 2 1*. Then (fTTl) and 
( |T8l) imply that > t> 3 and Wi > W3. Since u^iOj > 0, it follows that 



□ 



1M\S M 2 l l = vSw = av\W\ + av\w 2 — bv^w^ > av\Wi — bv^w^ > 0, 

as desired. 

We now calculate further: 



FE 



6 4 
6 4 
6 3 



{FEf 



F 2 E 



60 36 
60 36 
54 33 



18 12 
18 12 
18 11 



(FE) 2 - F 3 E 



54 36 
F 3 E= 54 36 
54 35 

6 ' 
6 
0-2 
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Applying Lemma [151 to S = (FE) 2 — F 3 E yields 

lM 1 (FE) 2 M 2 l t > lM 1 (F 3 E)M 2 l t 

for all Mi, M 2 G V, and hence, if F 3 E occurs in T{F 2 ■ ■ ■ T n , we can replace it with 
(FE) 2 without decreasing the norm of the product. Repeating this as many times as 
needed we can thus reduce the problem to products of the form (fT9l) with j\ £ {1, 2} 
for I = 1, . . . , k. 

It remains to eliminate factors F 2 E. Without loss of generality we may assume 
that the product T{F 2 ■ ■ ■ T n contains an even number of such factors. (Otherwise, 
we can left multiply T{T 2 • • ■ T n by F 2 E to create an even number of occurrences; 
the extra factor has no impact in the limit as n — > oo.) Note that two consecutive 
occurrences of F 2 E are separated by some power of FE. 

Lemma 16. For each k > 0, the matrix 

S k := {FE) k+3 - F 2 E(FE) k F 2 E 



satisfies the condition of Lemma {J 

Proof. Let G = FE. The characteristic polynomial of G is A 3 — 9A 2 — 6A, so G 
has three distinct real eigenvalues and p{G) = a. Furthermore, G satisfies its own 
characteristic equation: 

G 3 - 9G 2 -6G = 0. (20) 

The key to showing that Sk has the required form is the apparent coincidence that 
F 2 E is "almost" equal to 3FE. Precisely, F 2 E = 3FE + D, where D = [d*,] is the 
matrix whose only nonzero entry is 0^33 = 2. Multiplying (12"0"1) by G k and putting the 
results together we obtain, after some elementary algebra, 



S, = 6G k+1 - 3G k+1 D - 3DG k+1 - DG k D. 



(21) 



Writing 



(12"T|) becomes 



G k 



c*k (3 k 
Ik 5 k 



Sk 



Qa k+1 
6a k+1 
-64+1-44 



It remains to verify that 60:^+1 > Q5 k +i + 44- But this follows easily by induction 
since {a k } and {5 k } both satisfy the recursion x k +i = Qx k + Qx k -i for k>2, and the 
inequality clearly holds for k = 0, 1, 2. □ 
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In view of the last lemma, we can replace F 2 E(FE) k F 2 E with (FE) k+3 without 
decreasing the norm of the product T{T 2 ■ ■ ■ T n . Applying this repeatedly we eliminate 
all factors F 2 E two at a time. Therefore, the extremal case is T1T2 • • - T n — (FE) n ^ 2 
(assuming without loss of generality that n is even), and this shows that p({E, F}) = 
a/ p(FE) = y^a, proving the Proposition. □ 



We can now complete the proof of Theorem [3j By Lemma [13] and Proposition 
T4l (or rather, its proof), there is a constant C such that c n + a n < Ca n ^ 2 for all n. 
Using Lemma [TD] it follows that for each n, the level set Lf(y) can be covered by at 
most C(x n l 2 intervals of length 4 _n . Hence, 

dim// Lf(y) < lim — ^ — ; }- = — - — = d. 

- n _,oo _ log(4-») log 16 



giving the stated upper bound. 

3.3 Intersection with lines of integer slope 

We end this section with a proof of Corollary HJ It uses the following simple lemma: 

Lemma 17. For any f &T V and for any m e Z, the partial function f\ m \ has slope 
m on exactly one interval I\ m \,j with < j < 2' m L 

Proof. We prove the statement for m > 0; the case m < follows by symmetry. The 
statement is obvious for m = 0. Proceeding by induction, suppose f m has slope m on 
I mJ . Then f m+1 has slope m + lon I m +i,2j if u m>j = 1, or on 7 m+1)2i+1 if u mJ = -1. 
The uniqueness of j follows easily by induction as well. □ 

Proof of Corollary [^} For convenience we regard each / e T v as a 1-periodic function 
defined on all of R. However, we keep the convention that Graph(/) = {(x, f{x)) : 
< x < 1}. Fix m and b, and put £ := £ m ^- We first show that 

«~ P h(/)nq<{^ l f f H (22) 

Assume first that m > 0. Let 

m-l 

fa) ■= - E ^^(2^) + 2^ / ^' < ^ < 

fc=0 
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and define the linear mapping \1> : M 2 — > M 2 by y) := (2~ m x, —m2~ m x + 2~ m y). 
Then \1/ maps Graph( /) onto Graph^Lo^-mi), and it maps £ onto the horizontal line 
£' : y — b/2 m . Since \I/ is linear and invertible it is bi-Lipschitz, and hence 



Observe that g G 7^, and if / G 7^ then g E T c . Thus, ( J221) follows for the case m > 
from the bounds of Theorems Q] and [3j And for m < 0, it follows by applying the 
above argument to — /. 

We next show that the upper bounds are attained for each meZ. We do this for 
m > 0; the case m < follows by symmetry. First, let / G T c . By Lemma [T71 there is 
< j < 2 m such that f m has slope m on I m j. Let x = j/2 m be the left endpoint of 
I m j. The function := 2 m {j{2~ m x + x ) — f m {2r m x + x )) is in 7^, so its graph 
intersects some horizontal line I in a set of dimension 1/2. Let $ : M 2 — >• R 2 be the 
affine mapping that maps Graph(/i) onto Graph(/|,j m .). Then $ maps I onto a line 
with slope m, and since $ is bi-Lipschitz, it follows that 

dim H (Graph( /) n £') > dim H (Graph(/U m J n I') = dim H ( Graph (h) Hi) = 1/2. 

Finally, we build a function g a T v which attains the bound in (!22j) . Assume 
again that m > 0. Let h G 7^ and y £ K be such that dim H L h (y ) = d* v . Define 
g eT v hy 



Let $ be the affine mapping that maps Graph(/i) onto Graph (g |,[o,2- m ])- Then $ 
maps the line y = y onto a line £' of slope m, so that 

dim H (Graph(s() nf) > dim H (Graph(^L[o,2-"*]) H £') = dim^ L ft (y ) = cZ*. 



4 The random case 

In this section we prove the results for the random case. 
4.1 Dimension of the zero set 

Proof of Theorem® Observe first that, if f n does not take the value anywhere on 
an interval I n j, then / itself will not vanish anywhere in I n j. This holds since f n > 
on I n j implies f n > 2~ n on I n j, and \f — f n \ < 2~ n by Lemma [TU1 So for each level 



dim H (Graph(/) D £) = dim H ( Graph (#L[o )2 -™]) n £') < dim H L g {b/2 m ). 




< x < 1. 



This completes the proof. 



□ 
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n, we need only consider intervals I n j on which the graph of f n has at least one point 
on the x-axis. 

We calculate the almost-sure dimension of Lf(0) n [0, 1/2]. By symmetry, this will 
equal the almost-sure dimension of £/(0). The idea is to represent L/(0) fl [0, 1/2] as 
the attractor of a Mauldin- Williams random recursive construction; see [20] . When 
p = 1/2, the slope of /„ on J n = [0, 2~ n ] follows a simple symmetric random walk, 
and as such it returns to zero infinitely often with probability one. This must happen 
at even times. If it happens for the first time at time 2n, then f 2n = on / 2rii o, the 
slope of / 2n on the adjacent interval / 2rii i must be ±2, and / 2 „ cannot take the value 

anywhere else in [0,1/2]. Since Hausdorff dimension is independent of scale, we 
may assume without loss of generality that n = 1, and the slope of / 2 on J 2 .i is 2. 
We put J® := [0, 1/2], and call the graph of / 2 on [0, 1/2] our "basic figure". 

Next, the graph of / will evolve independently on the intervals [0,1/4] and 
[1/4,1/2]. After random waiting times T\,T% and T 3 , independent of each other, 
the basic figure will reappear at smaller scale (and possibly reflected in the x-axis, 
reflected left-to-right, or both) just to the right of 0, just to the left of 1/4 and just 
to the right of 1/4, respectively. Here T\ and T 2 have the distribution of 1 + n, and 
T3 has the distribution of r 2 , where for m G N, r m is the hitting time of the random 
walk of level m. For i = 1,2, 3, let Jj denote the projection of the ith copy of the 
basic figure onto the x-axis, so | Jj|/| J$\ = 2~ Ti . We note that outside J\ U J 2 U J3, / 
cannot take the value 0. 

This process continues. Each interval Jj generates, with probability one, three 
random nonoverlapping subintervals Jn, Jj 2 and above (or below) which the basic 
figure reappears for the second time in the construction of /, etc. This way, we obtain 
a ternary tree {J a : a G £*} of intervals, where S := {1,2,3} and S* := U^Lo^" - - 
The contraction ratios R aifi := | J a *i\/\Ja\, °" G E*, i = 1,2,3 are all independent. If 

1 G {1, 2}, Ro-^i has the distribution of 2~ ( - 1+T1 \ while has the distribution of 2~ T2 . 
Let K := fl^Li Uo-eS" ^ e the limit set. A by now familiar argument shows that 
K = L f (0) n [0, 1/2]. It follows from Theorem 1.1 of Mauldin and Williams [20] that 
with probability one, dim# K is the unique number s such that E(i?^ + R^ + = 1, 
in other words, the unique s such that 

2E (2- s ( 1+T1 )) + E (2- ST2 ) = 1. (23) 

Now let ipi(x) := E(x Ti ) denote the probability generating function of 7$. From stan- 
dard random walk theory, we have ^2 0*0 — {i ) i{x)) 2 , and V'i( a; ) = (1 — VI x 2 )/x. 
Putting r = 2 _s , fT23|) therefore becomes 

2r^i(r) + (Mr)) 2 = 1- 
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This equation has only one positive solution, given by r 2 = (V5- l)/2. Thus, 



logr _ log(r" 2 ) _ log ((1 + y/E)/2) _ 

— — «0j 



log 2 log 4 log 4 

and the proof is complete. □ 

Proof of Proposition^ The Gray Takagi function, defined by (j3J), satisfies the func- 
tional equation 

f(x)= U + |/(2*), 0<x<l/2 
\l -x- |/(2-2x), 1/2<z<1. 

Let x m := 1 - Ya=i A ~ (2i ~ 1) for m G Z +' and let x * '■= lim m-^oo^m- Applying flU} 
repeatedly it may be seen that for each m G N, 

ze[x m ,x m _i] =► f(x) = -4^ 2m -Vf(4 2m - 1 (x m ^ 1 -x)). (25) 

(The somewhat cumbersome algebraic details are omitted here, but (|25|) is best 
understood graphically.) Furthermore, f(x) > for < x < x*. As a result, Lf(0) 
consists of 0, x*, and an infinite sequence of nonoverlapping similar copies of Lf(0) 
itself, the mth copy being scaled by x m _i — x m = 4"( 2m ~ 1 ^ and reflected left-to-right. 
Thus, dim#L/(0) is the unique number s satisfying the Moran equation 

oo 

^ 4 -(2m-l)s = ^ 
m=l 

A routine calculation gives s = do. □ 

In the setting of Model 1, the theorem of Mauldin and Williams is not applicable 
because the random contraction vectors (i? CT *i , R a *2, Ra*3), cr G S* are neither in- 
dependent nor identically distributed. We must therefore work considerably harder 
just to obtain the inequalities of Theorem [51 

Proof of Theorem{5\(i). Assume without loss of generality that p > 1/2. Define 
random index sets 



n 

r + 

n 

rr 



T n := {j : < j < 2 , f 2n (x) = for some x G I 2n ,j}, 

{j eT n :f 2n >0 on I 2nJ and \ft n \ > on I° n j }, (26) 
{j er„:/ 2 „<0 on I 2n ,j and \f 2n \ > on I% nj }. 



Let K n := U Jg r n ^n,j an d K := n^Lo-^"- Then K = Lf(0). Let Q n denote the 
restriction of the graph of f 2n to K n , and for an interval / C [0, 1], let Q n [i denote 
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Figure 5: The two-step evolution of the Z-shape (top) and the cup-shape (bottom) 

the further restriction of Q n to K n D /. Suppose Y n contains three successive integers 
j,j + 1 and j + 2, and let J := I 2n j U hnj+i U ^2n,i+2- If t ne slopes of / 2n on the 
intervals hnj, hnj+i and hnj+2 are 2,0,2 or —2,0,-2, respectively, we call G n \j 
a Z -shape. If instead the slopes are —2, 0, 2, we call Q n [j a cup-shape. Let _/V n 
denote the total number of Z-shapes contained in Q n . We shall show that, given 
that N no > for some no, N n grows at an exponential rate with probability one. 

Assume iV no > 0. The top part of Figure [5] makes clear that if Q n [j is a Z-shape, 
then G n +ilj contains again a Z-shape, at 1/4 the scale of the original one. Hence, 
N n +i > N n . Moreover, if (cu2 n , 0J 2n +i) = (1, — 1), then Q n+ \[j contains in addition 
a cup-shape. As the bottom part of Figure shows, a cup-shape in Q n+ \ induces 
a pair of Z-shapes in Q n+ 2 if oo 2n +2 = — 1- Thus, a Z-shape in Q n induces three 
Z-shapes in Q n+2 if {oj 2n ,uj 2n+ i,uj 2n+2 ) = (1, —1, —1). By symmetry, the same is true 
if (uJ 2n ,u 2n+1 ,ui 2n+2 ) = (-1, 1, 1). As a result, P(iV n+ 2 > 3iV„) > pq 2 + qp 2 = pq. It 
now follows that for each m G N, A^ no+2m > N m N m , where N m is the product of m 
independent random variables Y\, . . . , Y m each having the distribution P(Yj = 3) = 
pq = 1 — P(Yi = 1). By the strong law of large numbers, 

logiV m 

> F(log Yi) = pqlogS a.s., 

m 

and since N n is nondecreasing, this implies that given iV no > 0, 

. , logA^„ pglog3 
hmmt > a.s. 

n->oo n 2 

Now fix e > 0. Given that N no > 0, there is with probability one an integer n\ 
such that 

logiV n pglog3 
n > n\ =>• > — e. (27) 

7?; _ 
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Choose 5 > so that 5 < 4~ ni , and let U be a covering of K = Lf(0) by intervals of 
length less than 5. By a standard argument, we may assume that U C {hnj '■ n > 
ni,0 < j < 2 2n }. Moreover, we may assume that if i2nj £ ^ j then j e T n . There 
certainly is a smallest n such that W contains some interval T 2 n,j\ fix this n for the 
remainder of the proof. If \f 2n \ > 4 on I 2n j, then r n+ i contains only one of the four 
integers 4j, . . . , 4j + 3 and a more efficient covering is obtained by replacing I 2n j 
with the corresponding subinterval (I 2n +2,Aj or I 2n +2^j+z)- Thus, we may assume 
that l/gj < 2 on I^j whenever I 2n j 6 £/. 

If ji and J2 are distinct members of T n such that /2 n = on hnji U ^2nj 2 > then 
L/(0) n/2n,ji and £/(0) PI hnjs are identical up to translation. Hence we may assume 
that U contains either all of the intervals I 2n j 011 which f 2n = (n still fixed), or 
none of these intervals. Similarly, if j\ and j 2 are distinct members of r+ such that 
\f 2n \ = 2 on I^ n ji U I 2n j 2 , then L/(0) n I 2n j 1 and £/(0) fl / 2 n,j 2 are identical up to 
translation and, possibly, a reflection. Hence we can assume that U contains either 
all of the intervals I 2n j with j G T+ on which f 2n has slope ±2, or none of these 
intervals. The analogous statement holds for j e r~. 

Now, since a Z-shape contains exactly one line segment of each of the three types 
considered in the last paragraph, it follows that the number of intervals I 2n j that lie 
in U is at least N n . Hence, 

J^^ 1 ^ - Nn ' A ~ US ( 28 ) 
leu 

for all s > 0. Taking 

f pq log 3 \ 
s = I — e I /log 4, 

we finally obtain from ff27|) and ( 128]) that Ylieu \^\ S — ^ Since e was arbitrary, this 
implies that 

dim^(0) > > 0, (29) 

almost surely given that N no > for some n . 

In order to complete the proof, we must show that 

P(jV„ = Vn) = 1 - ^ (30) 

and 

P(#L/(0) < oo|iV n = Vn) = 1. (31) 

Let <S n := Wo + • • • + w n -i- Then {S n } is a simple random walk with parameter p, 
and since p > 1/2, S n takes positive values infinitely often with probability one, but 
will visit the value —1 only with probability q/p. Observe now that S n is the slope 
of f n on J n = [0,2~ n ]. Suppose that for some integer n, S 2n -\ = — 1, S 2n = and 
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S 2n +i = 1- Then the slope of f 2n -i on / 2 n-i,o is — 1; the slope of f 2n is on I 2n ,o 
and -2 on I 2t1i1 ; and the slope of f 2n +i is 1 on J 2n +i,o and -1 on J 2n +i,i and I 2n+ i )2 . 
Now one sees that, regardless the value of u 2n +i, a Z-shape appears in Q n+ \ (part of 
the graph of / 2n + 2 ) somewhere on h n +i,i U hn+i,2- On the other hand, if S n > for 
every n, then f n > for each n and no Z-shape ever appears. Thus, 

P(iV n > for some n) = P(3n : S 2n -i = -1, = 0, S^+i = 1) = q/p, 

by the properties of the random walk. This gives (I3"U1) . 

Finally, note that the probability that S n > for all n and S n = for infinitely 
many n is zero. Thus, 

P(S n > for all but finitely many n\N n = Vn) = 1. 

But if S n > for all n and 5* n = for only finitely many n, then #T n is eventually 
constant. More precisely, there is no such that for every n > no, T n+ i contains 
exactly one integer from {4j, . . . , 4j + 3}, for each j e T„. This clearly implies that 
the limit set K = Lf(0) is finite. This proves ( 13TT) . and completes the proof of the 
theorem. □ 

Remark 18. The lower estimate ( 1291) for dimn Lf(0) can be improved by considering 
the ratio N n+ k/N n for values of k larger than 2, because there is a variety of ways 
for additional Z-shapes to appear in Gn+k, and this number increases exponentially 
as k increases. The details are cumbersome, however. Using a computer program 
with the value k = 9, the author has been able to establish, for the case p = 1/2, 
that dimff Lf(0) > .2548 a.s. This is still well below the upper bound of d ps .3471. 
The lower bounds appear to converge extremely slowly as k — > oo, and it is merely 
a guess that they converge to do. 

Proof of Theorem\B(ii) . Fix, for the time being, an integer k > 3. Recall the random 
index sets defined by ( 126]) . and define random variables 



' n 



= #{j eT n : f 2n = on I 2n ,j}, 

= #{j G T+ : |/^| = 2z on J 2 ° nj }, 1 < i < k, 

= #{j G r; : \& n \ = 2i on J° nJ }, 1 < i < k, 

= #{J G r+ : |/iJ > 2fc on J 2 ° nj }, 

= #0 G r; : |/^| > 2A; on J 2 ° nj }. 



The dynamics of these sequences of random variables depend on four cases regarding 
the signs (u 2n , w 2 „+i) , as follows. 
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Let z n := (cj, 
the random vectors z 1; . . 
four cases above all occur with probability 1/4, and hence we have 



,4 ), and let T n be the a-algebra generated by 
- 1, . . . , fc. Since p = 1/2, the 



Let a n l) := Un' + In 



(<) , 7 W 



E(c 



n+l 



(0 

n+l 



E(<7, 



E(a (fc) 



n+l 



l^n) 
l^n) 



n i ^ n 5 

4 ™ 



2c r , 



4°" ' 



1 



(32) 



2 < i < fc, 



Put x n := (c„, <r n } , . . . , a n ') . By 

E(x n+1 |J r n ) < A fc x n , 

where is the (fc + 1) x (fc + 1) tridiagonal matrix 

11/4 

2 3/4 1/4 '•• 
1/4 1/2 1/4 ••• 



Ai 







1/4 1/2 1/4 
1/4 1 



(33) 
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Let pk denote the spectral radius of Ak- Since Ak is nonnegative pk is an eigenvalue 
of Ak, and since Ak is irreducible, the Perron- Frobenius theorem guarantees the 
existence of a positive left eigenvector of Ak corresponding to pk- It follows by 
( 1331) that the process 

x n ■= Pfe n VjfcX n , neZ+ (34) 

is a positive supermartingale, which by the Martingale Convergence Theorem con- 
verges almost surely to a finite nonnegative limit X^. Let 5 be the smallest entry of 
the vector Vfc. Then 5 > 0, and for a given e > 0, (1341 implies that for all sufficiently 
large n, 



C + ^a«<(l + e)r 1 X 00 ^. 



8=1 



Thus, the number of intervals l2 n ,j needed to cover Lf(0) grows at most at rate pjj. 
Consequently, 

a- t in\ ^ lo SPfc 
dmitf LfiO) < a.s. 

log 4 

The proof will be complete if we can show that 

liminf p k < 1 (35) 

k— >oo Z 

Let Ak be the fcxl; matrix obtained by deleting the last row and last column of Ak- 



For completeness, we define A 1 = [1] and A 2 



1 1/4 

2 3/4 



Let 



&(A) := det(A fc - A/ fc+ i), Cfc(A) := det(A fc - XI k ) 
be the characteristic polynomials of Ak and Ak, respectively. Then 

&(A) = (l-A)C*(A)-ic*-i(A), (36) 

lb 

and Cfc satisfies the recursion 

C*(A) = Q " A ) Cfe-i(A) - ^C fc - 2 (A), fc > 3, (37) 

with initial conditions 

Ci(A) = l-A, C 2 (A) = A 2 -^A + i (38) 

For fixed real A > 1, the solution of the system (|37|) . (|38|) can be written as 

G(A) = C 1 (A) 7 f- 1 (A) + C 2 (A) 72 fc - 1 (A), (39) 



28 



where 

7i(A) :=iQ-A- v/A^a) , 72 (A) := \ Q - A + v/A^V) , (40) 

and 

r m _ Ci(A) 72 (A)-C 2 (A) / n _ C 2 (A)-Ci(A)7i(A) , . 

72(A) -7l(A) 72(A) -7l(A) 

Substituting f )39|) into ( 136]) and some rearranging finally gives 

6(A) = C 1 (A) 7l fc - 2 (A) + C 2 (A) 7 r 2 (A), fc > 3, (42) 
for real A > 1, where 

a(A):= ((l-A) 7i (A)-^C i (A), i = 1, 2. (43) 

Now observe that A > 1 implies (72(A)! < 1/4 and 71(A) < 0, with |7i(A)| > 1 if 
and only if A > 25/16. If liminf p*. < 25/16 we are done, since 25/16 = 1.5625 < 
(1 + y/o)/2. So assume liminf pf, > 25/16; then we can find p > such that 
Pk > 25/16 + 77 for all sufficiently large k, and hence |7i(pfc)| > 1 + e for all sufficiently 
large k, for some e > depending on 77. Since £fc(Pfc) = for all k, this, together 
with ( ]42l . implies that Ci(pfc) must tend to zero as k — > 00. Hence, by ([431 . either 
Ci(Pfe) -)• or (1 - Pk)li(Pk) 1/16. It is easy to see that the latter is impossible, 
and so p* := liminf p^. must be a solution of Ci(p*) = 0, as Ci(A) is continuous in A. 
Routine algebra using (1381) . (j4"0l and (jUJ) shows that the only solution of Ci(p*) = 
is p* = (1 + v / 5)/2. This proves ( 1351) . completing the proof of the theorem. □ 



4.2 Dimension of the maximum set 

Proof of Theorem^ Slightly abusing notation, put M n := max{/ n (x) : < x < 1} 
and M n := {x G [0, 1] : f n (x) = M n }. Note, since <f>(x) + (l/2)0(2x) < 1/2, that 

n_1 1 /IV 

i=0 ^ ' 

Define random index sets 

A^:=|j:0<i<2 2n ,/ 2n = ^^ on/ 2n ,A, n E Z + , 

and put X n := #A*. Then {A 7 "*} is a Galton- Watson process with initial value 
Xl = 1 and offspring distribution 7r = (7r , 7Ti, 7r 2 ) = (1 — 2p 2 g — p 3 , 2p 2 g,p 3 ), where 
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q = 1 — p. The offspring distribution has mean p := Y2i=o'i 1T i = ^P 2 and probability 
generating function h(t) := ^Lo 71 "*^ = 0- ~~ ^P 2< i — P 3 ) + 2p 2 qt + pH 2 . Let p : = 
P(X* — > 0). According to the basic theory of the GW process (e.g. [9j), p = 1 when 
p < 1, and in that case, Mf < 2/3 with probability one. Assume from now on that 
p, > 1; that is, p > l/y/2. Then p is the smallest positive number satisfying t = h(t), 
so that p = (1 — 2p 2 g — p 3 )/p 3 , and 

2d 2 — 1 

P(M, = 2/3) = P(X* >0Vn) = l-p= ^ . 
This establishes part (i) of the theorem. Next, define the random set 

oo 

■■ n u w 

n=0j£Ai 

Then = if and only if X* — )■ 0, and given that F\ ^ 0, dim# F x = log/// log4 
a.s. by Theorem 1.1 of Mauldin and Williams [20J. 

Put t = 0. Proceeding inductively, suppose processes {X*}, ■ ■ ■ , {X*} and ran- 
dom variables t^, ... , t^-i have been defined and that X r n — » for r = 1, . . . , k. Let 

:= min{n : X% = 0}, and define 

Tfc := inf{n > Tk-i + A& : M.2n contains an interval}. 

Since p > 1/2, T k is finite almost surely. (Put m = r fe _ 1 + N k , and let x ■= j/2 2m 
be a point of maximum of fim- The slope of fim+n directly to the right of xo starts 
with a nonpositive value and follows (as a function of n G Z + ) a simple random walk 
with parameter p, so it will eventually reach 0.) Note that 



n-l Wi \ T-fc+« 

M 2(rfe+n) <M 2r ,+J]-( I ] . „eS H 

i=0 



Define the random index sets 



f n_1 1 /l\ Tk+i 

K +1 ■= <3 : < j < 2 2 ^\f 2(Tk+n) = M 2Th +E2UJ ° n 



for n G Z+, and put X^ +1 := #A£ +1 . By definition of r fc , A k +l ^ and so X fc+1 > 1. 
Now {X^' +1 } is again a Galton- Watson process with offspring distribution 7r, and it 
depends on the preceding processes {X*},...,{X*} only through the value of X fc+1 
Thus, we have 

P(X n fc+1 > \/n\X l n ^ 0, . . . , X k n ^ 0) 

> P(X n fc+1 > \/n\Xl^ 0, . . . , X k n 0, X fc+1 = 1) = 1 - p. 
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Define the random set 

oo 

Fk+i == P| [J h(Tk+n),j- 

Then F k+ i = if and only if X^ +1 — > as n — > oo, and given that X^ +1 > for all 
n, dim H F k+1 = log /x/ log 4 a.s. 

With probability one, there will eventually be a k G N such that X% > for all 
n, and for that k, we have F k = M./, and 

Part (ii) of the theorem now follows. □ 

5 Open problems 

There are many natural questions left to answer. A few are listed here. 

Problem 1. Does there exist / G T v such that dim^ Lf(y) = for every 

Problem 2. Is it true for all / G T v that Lf(y) is finite for Lebesgue-almost every 
yl Is this true for the Gray Takagi function? 

Problem 3. (Random case, Model 1) Prove or disprove that dim# Lf(0) = do a.s. 
when p — 1/2. 

Problem 4. (Random case, Model 2) What is the best (smallest) bound Si such 
that, for each y e R, P(dim^ Lf(y) < s\) = 1? More strongly, what is the smallest 
s 2 such that P(dim^ Lf(y) < s 2 Wy G E) = 1? (Obviously, s a < s 2 < d*.) In the 
case of Model 1, only the first question is of interest, as s 2 — 1/2 in view of Theorem 

m 

Problem 5. (Random case, Model 2) Prove or disprove that dim^A4/ = a.s. 
when p < 1 / "\/2, and that A4f is finite a.s. when p < 1/2. What can one say about 
the distribution of M/ in these cases? 
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